Bell inequalities for continuous- variable correlations 
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We derive a new class of correlation Bell-type inequalities. The inequalities are valid for any 
number of outcomes of two observables per each of n parties, including continuous and unbounded 
observables. We show that there are no first-moment correlation Bell inequalities for that scenario, 
but such inequalities can be found if one considers at least second moments. The derivation stems 
from a simple variance inequality by setting local commutators to zero. We show that above a 
constant detector efficiency threshold, the continuous variable Bell violation can survive even in the 
macroscopic limit of large n.This method can be used to derive other well-known Bell inequalities, 
shedding new light on the importance of non-commutativity for violations of local realism. 



Einstein, Podolsky and Rosen (EPR), in their famous 
1935 paper [3], demonstrated the incompatibility be- 
tween the premises of local realism and the completeness 
of quantum mechanics. The original EPR paper used 
continuous position and momentum variables, and relied 
on their commutation relations, via the corresponding 
uncertainty principle. Bohm Q introduced, in 1951, his 
version of the EPR paradox with spin observables. This 
was the version that was used by Bell 0] to prove his fa- 
mous theorem showing that quantum mechanics predicts 
results which can rule out the whole class of local hid- 
den variable (LHV) theories. It is hard to overemphasize 
the importance of this result, which has even been called 
"the most profound discovery of science" However, 
the original Bell inequality, and all of its generalizations, 
are directly applicable only to the case of discrete ob- 
servables. The main purpose of this letter is to close the 
circle and derive a class of Bell-type inequalities appli- 
cable to continuous variables correlations, together with 
multipartite generalizations. 

We derive a class of inequalities for local realism that 
directly use correlations of measurements, with no re- 
striction to spin measurements or discrete binning. The 
new inequalities are remarkably simple. They place no 
restriction on the number of possible outcomes, and the 
contrast between the classical and quantum bounds in- 
volves commutation relations in a central way. They 
must be satisfied by any observations in an LHV theory, 
whether having discrete, continuous or unbounded out- 
comes. We can immediately rederive previously known 
Bell- type inequalities, obtaining at the same time their 
quantum-mechanical bounds by considering the non- 
commutativity of the observables involved. We also dis- 
play quantum states that directly violate the new in- 
equalities for continuous, unbounded measurements, even 
in the macroscopic, large n limit @, USB- We show 
that the new Bell violations survive the effects of finite 
generation and detection efficiency. This is very surpris- 
ing, in view of the many examples in which decoherence 
rapidly destroys macroscopic superpositions [§]. 

Apart from this intrinsic interest, these inequalities are 
relevant to an important scientific problem. No experi- 
ment has yet produced a Bell inequality violation without 



introducing either locality or detection loopholes. One 
path towards this goal is to use continuous-variables (CV) 
and efficient homodyne detection, which allows much 
higher detection efficiency than is feasible with discrete 
spin or photo-detection measurements. A number of 
loop-hole free proposals exist in the literature, but they 
all use Bell [lfl, ELI H El Q or Hardy [H] inequalities 
with a dichotomic binning of the results (which usually 
lead to small violations), or else a parity or pseudo-spin 
approach[l(| E3, EI| which cannot be realized with ef- 
ficient homodyne detection. Are there Bell inequalities 
which can be derived without the assumption of finite 
number of outcomes and therefore are directly applica- 
ble to CV - with no need to bin the results? 

For n parties, m measurements per party and o out- 
comes, it is well-known that the set of correlations al- 
lowed by LHV theories can be represented as a con- 
vex polytope, a multi-dimensional geometrical structure 
formed by all convex combinations (linear combinations 
where the coefficients are probabilities, i.e., they are non- 
negative and sum to one) of a finite number of vertices. 
The vertices of this polytope are the classical pure states 
— the states with well-defined values for all variables 
0, EH The tight Bell inequalities are associated 

with the linear facets of the polytope. It is a computa- 
tionally hard problem to list all Bell inequalities for given 
(n, m, o), and full numerical characterizations have been 
accomplished only for small values of those parameters. 

However, no class of Bell inequalities has previously 
been derived without any reference to the number of 
outcomes or to their bound. Any real experiment will 
always yield a finite number of outcomes; but are there 
constraints imposed by LHV theories that are indepen- 
dent of any particular discretization, and can be explic- 
itly written even in the limit o — > oo? Our answer is yes; 
and the derivation is much more straightforward than 
in the case of the usual Bell-type inequalities, which are 
restricted to a particular set of outputs. 

We will focus on the correlation functions of observ- 
ables for n sites or observers, each equipped with m pos- 
sible apparatus settings to make their causally separated 
measurements. We consider any real, complex or vector 
function F (X, Y, Z, . . .) of local observations Xi, Yi, Zi 
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at each site i, which in an LHV theory are all functions 
of hidden variables A. In a real experiment the different 
terms in F may not all be measurable at once, because 
they may involve different choices of incompatible observ- 
ables. The assumption of locality enters the reasoning by 
requiring that the local choice of observable does not af- 
fect the correlations between variables at different sites, 
and therefore that the averages are taken over the same 
hidden variable ensemble P (A) for all terms. We intro- 
duce averages over the LHV ensemble (there's no loss of 
generality in considering deterministic LHVs (2ll|). 

(F)= J P(X)F(X(X),Y(X),Z(X),...)dX . (1) 

Our LHV inequality uses the simple result that any 
function of random variables has a non-negative variance, 

|(F)| 2 <(|F| 2 ). (2) 

We can also give a bound (|-F| 2 ) < (\F\ 2 ) sup , where the 
subscript denotes the supremum (least upper bound), in 
which products of incompatible observables are replaced 
by their maximum achievable values. This is necessary 
since if we are not able to measure both Xi and Yi si- 
multaneously, a general LHV model could predict any 
achievable correlation [22l |. 

The same variance inequality applies to the corre- 
sponding Hermitian operator F in quantum mechanics. 
While the observables at different sites commute — they 
can be simultaneously measured — those at the same 
site do not, so operator ordering must be included. This 
enables us to see how quantum theory can violate the 
variance bound for an LHV. 

As an example, we will apply this variance inequality to 
a well-known case. Consider two dichotomic observables 
Xi, Yi per site i, the outcomes of which are ±1. We define 
F\ = X\ , F{ = Y\ , and then inductively construct [23| : 

Fn = 2 (F n -i + F'i-i) x n + ^{Fn-i - F' n _i)Y n , (3) 

where F n can be obtained from F n by the exchange 
Xi < — ► Yi. In calculating F% we'll keep track of the lo- 
cal commutators just to make the contrast with quantum 
mechanics clearer. For real variables X, Y, the commu- 
tator is defined in the same way as for the corresponding 
operators, i.e., [X, Y] = XY—YX. The anti-commutator 
is defined by [X, Y]+ = XY + YX. Then 

F n = 2 {C*n-1 + Fn-l)(Xn + Y n) 

+ [F n -i,F n _ 1 ] + (X n - Y n ) + (F^_ 1 - Fl l _ l )[X n ,Y n ]+ 

-[F n _ x ,F^_ t }[X n ,Y n ]}. (4) 

Since X n = Y% = 1, we can show that F 2 = F n 2 and 

Fl = Fn-l - \ [Fn-1, K-i] [Xn, Y n ] . (5) 



In a LHV theory, the term which involves commuta- 
tors will be zero since [X(X),Y(X)] = X{X)Y(X) - 
Y(X)X(X) = 0. Hence by induction F% = f} = 1 and 
the variance inequality ((2j becomes: —1 < (F n ) < 1. 
This is the Mermin-Ardehali-Belinskii-Klyshko (MABK) 
0, [H, Bell inequality, which reduces to the well- 
known Bell-CHSH [2a| inequality for n = 2. 

We can now calculate the quantum mechanical bound 
by writing the variance inequality J2]) and substituting 
the functions in <[5j) by their corresponding operators 

(F n ) 2 Q < (Fl) Q ^(Fl_ 1 -\[F n ^ 1 ,F n _ 1 ][X n X]) Q 

< \\FLl\\ + {\\[h-lJLl}\\\\[XnX}\\, (6) 

where the norm \\A\\ denotes the modulus of the maxi- 
mum value of (A)q over all quantum states. The norm 
of the second commutator has the bound || [X n , Y n ]\\ < 
2. It's easy to show that [F n ,F n ] = F^lXnXn] + 
[Jn-i.^-i] and therefore \\[F n ,F^\\ < 2\\F? l _ 1 \\ + 
||[^ra-i,i^_x]|| • Solving the recursion relation by not- 
ing that \\F^\\ = i||[A 1 ,Y 1 ]|| = 1 we finally arrive at 
the bound (F n )Q < 2 n ~ 1 . This can be attained with the 
generalized GHZ states j23|, which therefore violate ([2]). 

Inspired by those results, we now demonstrate an LHV 
inequality that is directly applicable to unbounded con- 
tinuous variables, in particular field quadrature opera- 
tors. The choice of the function F n in ([3]) is not optimal 
though, since the variance in general involves incompat- 
ible operator products that have no upper bound. 

To overcome this problem, consider a complex function 
C„ of the local real observables {Xk, Yk} defined as: 

n 

C n = X n + iY n = 11 (X k +iY k ), (7) 
fc=i 

so that the modulus square only involves compatible op- 
erator products, i.e. \C n \ 2 = T[k=i( X k + Y k ) ■ Applying 
the variance inequality to both X n and Y n , we find that: 



(X n f + (Y n f<(H(Xl + Y 2 )) (8) 
fe=i 

This is our main result. Given the assumption of local 
hidden variables, this inequality must be satisfied for any 
set of observables Xk, Yk, regardless of their spectrum. 

The fact that we have neglected the commutators in 
deriving ^ hints that quantum mechanics might predict 
a violation. We define quadrature operators 

X k = a k e~ Wk + a\e^ 

Y k = a fe e" 4(efc+Sfc7r/2) +4e i(9fc+W2) , (9) 

where a k , d\ are the boson annihilation and creation op- 
erators at site k and Sk € { — 1, 1}. 
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We now define the operator Z k = X k + iY k and 
note that it follows that C n — TTk=i %k- The defini- 
tion of Yfc allows for the choice of the relative phase 
with respect to X k to be ±ir/2. Depending on s k , 



for each k either Zk 



2a k e 



or Z k 



2ale ie >>. 



Denoting A k (l) = a k and Afc(— 1) = aj,, the term 
in the LHS of iJH) in quantum mechanics is then 
IOL^qI 2 = \2 n e^^{Yl k A k {s k )) Q \\ The RHS 
becomes (Ilfe^i^Q'fcflfc + 2))q regardless of the phase 
choices. To violate ((Sj) we must therefore find a state 
that satisfies 



> 



n« 



a\a k 



(10) 



which is surprisingly insensitive to relative phases be- 
tween the quadrature measurements at different sites. 

This violation of a continuous variable Bell inequality 
can be realized within quantum mechanics. Consider an 
even number of sites, choosing s k = 1 for the first half of 
them and s k = — 1 for the remaining. To maximize the 
LHS we need a superposition of terms which are coupled 
by that product of annihilation/creation operators. One 
choice is a state of type 



|* s > = co |0, . . . , 0, 1, . . . , 1) + ci 1 1, . . . , 1, 0, . . . , 0> , 

(11) 

where in the first term the first n/2 modes are occu- 
pied by zero photons and the remaining by 1; conversely 
for the second term. With that choice of state the 
LHS of lflO|) becomes |co| 2 |ci| 2 , which is maximized by 
|c | 2 = |ci| 2 = f. The RHS is (§)*(§)* independently 
of the amplitudes cq, c%. Dividing the LHS by the RHS, 

inequality (flOf becomes \ (| ) 2 < 1, which is violated 
for n > 10, and the violation grows exponentially with 
the number of sites. 

While setting up the homodyne detectors necessary for 
this observation is challenging, the complexity of this task 
scales linearly with the number of modes. A more strin- 
gent constraint is most likely in the state preparation, but 
we can relate state i|ll[l to a class of states of great experi- 
mental interest. They can be achieved from a generalized 
GHZ state of n/2 photons, -^(|i/)®f + |F)®t ) — where 

| if) and \V) respectively represent single-particle states 
of horizontal and vertical polarization — by splitting each 
mode with a polarizing beam splitter. Therefore viola- 
tion of ([E]) can be observed in the ideal case with a 5-qubit 
photon polarisation GHZ state and homodyne detection. 

An interesting question is the effect of decoherence, 
both from state preparation error [22J and detector in- 
efficiency. The usual Bell-CHSH violations have an ef- 
ficiency threshold [28| of 83%. This has not yet been 
achieved for single-photon counting. Homodyne detec- 
tion is remarkably efficient by comparison, with up to 
99% efficiencies being reported. However, the effect of de- 
tector efficiency is easily included by assuming that each 




Figure 1: Minimum state preparation fidelity e m in for ideal 
detectors (solid line), and minimum detection efficiency rjmin 
for ideal state preparation (dashed line) required for violation 
of ([8]) as a function of the number of modes. The asymptotic 
value of r/min is indicated by the dash-dotted line. 



detected photon mode is preceded by a beamsplitter with 
intensity transmission n < 1. This changes both the LHS 
and RHS, so that the inequality becomes 2n+l — 4 2 ^", 
giving a threshold e fficiency re quirement of n > rj m i n , 
where r] min = (1 + %/l + 4 1 - 2 /«)/4 1 " 2 /". 

This reduces at large n to an asymptotic value of 
Voo = 0.80902. Unexpectedly, the Bell violation (which 
signifies a quantum superposition) is less sensitive to de- 
tector inefficiency in the macroscopic, large n limit. The 
minimum detector efficiency r\ n at finite n is plotted in 
Fig. 1, together with the minimum state preparation fi- 
delity t m in in the case of ideal detectors, where we model 
the density matrix as p = e\$>s){^s\ + (1 — e )f- 

We will finally prove that there are no LHV inequalities 
possible if one considers only the first-moment correla- 
tions between continuous variables in different sites. We 
will show this explicitly for the simplest case and indicate 
how to generalize to arbitrary numbers of parties and set- 
tings. Consider first n = 2 parties, Alice and Bob, each 
of which can choose between m = 2 observables: X a ,Y a 
for Alice and Xb,Yb for Bob. Each measurement yields 
an outcome in the real numbers. The first-moment corre- 
lation functions for each of the 4 possible configurations 
are just the averages (X a X b ), (X a Yb), (Y a X b ), (Y a Y b ). 
Given those 4 experimental outcomes, can we find a lo- 
cal hidden variable model which reproduces them? 

We construct an explicit example. Consider a hidden- 
variable state S where the hidden variables are the mea- 
sured values X, Y, in an equal mixture of four classical 
pure states S k — (X ai Y a , X b ,Y b ) k defined by 



51 = 2(1,0, (X a X b ),0) 

5 2 = 2(l,Q,0,{X a Y b }) 

5 3 = 2(0,l,{Y a X b },0) 

5 4 = 2 (0,1,0, {Y a Y b }). 



(12) 
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Each of the states Sk assigns a nonzero value to only 
one of the 4 correlation functions. Since the probability 
of each of the states in the equal mixture is 1/4, we have 
for example (X a X b ) s = ^.(ttjs, = (X a X b ). 

Satisfying the two-site correlations using the state S 
defined by (fl2|) leaves us with uncontrolled values for the 
single-site correlations, for instance (Xb)s = ^((X a Xb) + 
(Y a Xb)). One might object to the fact that this is not 
equal to (Xb) in general. However, we may correct these 
lower order correlations by adding four more states (S 5 
to Ss) and changing the prefactors multiplying Si to 

54 to compensate for their reduced weight in the equal 
mixture. Crucially, adding these extra states to S in 
this manner does not modify the values of correlations 
such as (X a Xb)- As an example, we exhibit the state 

5 5 = 8 (0, 0, (X b ) - ((X a Xb) + (Y a X b ))/VS, 0), which cor- 
rects the single expectation value (Xb)s to (Xb). 

The proof generalizes easily to arbitrary n and m. In 
that case, there are m n possible combinations of measure- 
ments which yield n-site correlations. Denoting the j th 
observable at site i by Xf, each combination is specified 
by a sequence of indices (ji,.]2, ■ ■ ■ ,jn)- For each combi- 
nation of measurements, we define a hidden variable state 
which assigns nonzero values only to the variables which 
appear in the associated correlation function (11™= 1 -^T)- 
In analogy to the example above, we can always choose 
the values of the hidden variables associated to X 3 ' such 
that their product is equal to m n (11™= 1 X 3z ) . Since all 
other m n — 1 states defined in this way will give a value 
of zero to this particular correlation function, and given 
that the probability associated with each of those states 
is l/m n , we reproduce all correlations as desired. As 



indicated in the example, additional first moment corre- 
lations involving less than n sites can be included in the 
LHV model by adding additional states to S in a way 
which doesn't affect the n-site correlations. Thus, any 
possible observation of first moment correlations may be 
explained using a LHV model, and hence these correla- 
tions alone cannot violate any Bell inequality. In other 
words, the minimum requirement for a correlation Bell 
inequality with continuous, unbounded variables, is to 
use not just the first but also the second moments at 
each site. 

In conclusion, we have derived a new class of Bell-type 
inequalities valid for continuous and unbounded experi- 
mental outcomes. We have shown that the same proce- 
dure allows one to derive the MABK class of Bell inequal- 
ities and their corresponding quantum bounds. That 
derivation makes it explicit that non-zero commutators 
— associated with the incompatibility of the local observ- 
ables — are the essential ingredient responsible for the 
discrepancy between quantum mechanics and local hid- 
den variable theories. The new Bell- type inequality de- 
rived here can be directly applied to continuous variables 
without the need for a specific binning of the measure- 
ment outcomes. Surprisingly, quantum mechanics pre- 
dicts exponentially increasing violations of the inequality 
for macroscopically large numbers of sites, even including 
realistic decoherence effects like inefficient state prepara- 
tion, and a detector loss at every site. 
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